Fermionic spin excitations in two and three-dimensional antiferromagnets 
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Spin excitations in an ordered Heisenberg magnet are magnons — bosons with spin 1. That may 
change when frustration and quantum fluctuations suppress the order and restore the spin-rotation 
symmetry. We show that spin excitations in the S = 1/2 Heisenberg antiferromagnet on kagome are 
spinous — fermions with spin 1/2. In the ground state the system can be described as a collection 
of small, heavy pairs of spinous with spin 0. A magnetic excitation of lowest energy amounts to 
breaking up a pair into two spinous at a cost of 0.06J. 



The Heisenberg antiferromagnet, described by the ex- 
change Hamiltonian 



(1) 



where J > is the exchange coupling and (ij) denotes 
a pair of neighboring sites, is a simple model of an with 
realistic prototypes. A bipartite antiferromagnet in three 
dimensions exhibits long-range spin order that breaks the 
global spin SU(2) symmetry. Low-energy excitations of 
such a magnet are spin waves, whose quantization yields 
magnons — bosons with spin 1 [I]. In one dimension, spin 
order is disrupted by quantum fluctuations and the SU(2) 
symmetry is restored. For spins of length S = 1/2, the 
excitations are spinous — quasiparticles with spin 1/2 [2]. 
A long-standing question is the existence of similarly un- 
usual excitations in higher-dimensional magnets where 
the symmetry of the ground state is restored by a combi- 
nation of strong quantum fluctuations (small S) and ge- 
ometrical frustration (a non-bipartite lattice). In 1973, 
Anderson proposed a resonating valence-bond state for 
spins 1/2 on the triangular lattice [3]. In this state, spins 
form singlet bonds, or quantum dimers [4], with their 
neighbors. Although the ground state on the triangular 
lattice turned out to be ordered, kagome and hyperk- 
agome lattices have emerged as likely candidates for ex- 
otic quantum physics [•")]. 

The exchange energy (1) on kagome [Fig. 1(a)] is min- 
imized when the total spin of every triangle attains its 
lowest allowed value [8]. For S = 1/2, that can be 
achieved by putting two spins of a triangle in a singlet 
state; the third spin is free to form a singlet on the adja- 
cent triangle. If a dimer could be placed on every trian- 
gle, we would construct a ground state with frozen sin- 
glets, a valence-bond solid. That does not work: one in 
four triangles on kagome lack a dimer. Quantum fluctu- 
ations induced by defect triangles make the ground state 
a nontrivial superposition of valence-bond states. Some 
recent theoretical works [9, 10] lend support to a valence- 
bond crystal proposed by Marston and Zeng [11], while 
others are consistent with a valence-bond liquid [12]. Ex- 
act diagonalization studies [lo] reveal a large number of 
low-lying singlet states, presumably associated with dif- 
ferent dimer configurations. Spin-1 excitations appear to 




FIG. 1: (a) Kagome lattice and the major players: quantum 
dimers (thick bonds), a defect triangle containing no dimer 
(red), and an antikink spinon (blue arrow), (b) The map- 
ping to a compact U(l) gauge theory on a honeycomb lattice. 
Electric flux on links has strength ±1 and is depicted as ar- 
rows in the manner of Elser and Zeng [(i, 7]. The inset shows 
the U(l) charges of the lattice sites and quasiparticles. 



have an energy gap of 0.05J to O.IOJ [12, 13, 14]. 

We propose a simple physical picture of the = 1/2 
kagome antiferromagnet, in which the system is viewed as 
an ensemble of spinous, fermions with spin 1/2. The anti- 
symmetry of the many-body wavefunction can be traced 
to quantum interference. As two spinous are exchanged, 
they drag around pairs of spins entangled in 5 = states. 
In the end, the background singlet pairs return to the 
same positions. However, an odd number of them reverse 
their orientation, thus altering the sign of the many-body 
wavefunction. 

The defect triangles turn out to be small, heavy pairs 
of spinous bound by exchange-mediated attraction. The 
binding energy is Ei, ~ 0.06J. Spin-0 excitations are 
associated with the motion of pairs; their heavy mass 
is reflected in the large density of singlet states at low 
energies. Spin-1 excitations correspond to breaking up a 
bound pair into (nearly) free spinous with parallel spins. 
Because spinon pairs retain their individual character on 
kagome, the spin gap is determined primarily by their 
binding energy, A w Eb = 0.06 J. This number is in line 
with the previous estimates of the spin gap [12, 13, 14]. 

This picture is based on a study of two toy models on 
lattices that share with kagome the triangular motif: the 
A chain [1-5, IG] and the Husimi cactus, a tree-like variant 
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FIG. 2: Ground states L (a) and R (b) of the A chain. A 
local triplet excitation (c) decays into a localized kink and 
a mobile antikink (d-e). Arrows on the dimers illustrate the 
pivoting rule (see text). 



of kagome [6, 17]. Both systems have dimerized ground 
states with exactly one singlet bond on every triangle. 
The A chain provides information about the properties 
of individual spinons, while the Husimi cactus sheds light 
on their quantum statistics and interactions. 

The A chain has two ground states, one with singlets 
on the left bonds of the triangles (L), the other with sin- 
glets on the right bonds (i?), Fig. 2(a-b). A local spin-1 
excitation decays into two spinons — quasiparticles with 
spin 1/2, Fig. 2(c-e). They serve as domain walls and 
come in two flavors [15]: kinks interpolate between state 
L on the left and state R on the right, antikinks do the 
opposite. A kink is fully localized and has zero energy 
cost: every triangle in its vicinity has a total spin 1/2 
and is thus in a state of lowest energy, which is of course 
a stationary state. One might think that zero excita- 
tion energy would lead to a proliferation of kinks, but 
that does not happen: kinks can only be created in pairs 
with antikinks, whose excitation energy is positive. That 
hints at the existence of a conserved charge, which will 
be defined later. 

Antikinks are mobile; their dynamics is described ap- 
proximately as hopping to the nearest triangle [15]: 
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where n is the triangle containing the antikink. A Fourier 
transform yields the quasiparticlc dispersion e(/c) = 
5 J/4 — J cos fc ~ J/4-1- A;^/ (2m) in the limit of low lattice 
momenta; the antikink mass m — 1/ J. A local spin-1 ex- 
citation results in the creation of a kink and an antikink 
with a minimum energy of J/4. Terms omitted in Eq. (2) 
create additional excitations by promoting the two sin- 
glets adjacent to the antikink to triplets; alternatively, 
they can be viewed as creation of an extra kink-antikink 
pair next to the antikink. These virtual excitations renor- 
malize the antikink mass to m w 1.16/J and lower the 
bottom of the antikink band to A = 0.219J [is] in agree- 
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FIG. 3: A dimerized state on the Husimi cactus. A local 
triplet excitation (a) decays into a localized kink and a mobile 
antikink (b-d). Displaced valence bonds are shown in blue 
color. Shaded triangles in (d) mark the allowed path of the 
antikink. 



ment with exact diagonalization [19]. Such perturbations 
appear to be harmless and will be ignored. 

We note that the sign of the hopping terms in Eq. (2) 
depends on the sign convention for quantum dimers [2i)]: 
the S = state of two spins is antisymmetric under ex- 
change: |(z,j)> = (I Ulj)-\ Uh))/V2 = -\U,t)). When 
necessary, we will depict a dimer state \ as an arrow 
pointing from i to j. The reversal of the arrow amounts 
to multiplying the state wavefunction by —1. A negative 
hopping amplitude is enforced by the pivoting rule: when 
an antikink hops from one triangle to another, a quan- 
tum dimer moves in the opposite direction by pivoting 
on the site shared by the triangles. Fig. 2(d-e). 

A ground state of the cactus has a quantum dimer on 
every triangle [17]. Like on the A chain, a local spin-1 ex- 
citation decays into a localized kink and a mobile antikink 
propagating along a one- dimensional path (Fig. 3). The 
kink has energy 0, while the antikink propagates along 
the allowed path with a minimum energy of J/4 and a 
mass TO = 1/J. (Again, we neglect virtual excitations in 
the vicinity of the antikink.) 

Next we examine a single defect triangle on the cac- 
tus. Fig. 4(a). The Hamiltonian connects this state to a 
state with a singlet bond of a longer range [G], Fig. 4(b). 
The longer-range bond can be viewed as two antikinks 
with total spin 0. By jumping to a neighboring tri- 
angle, the antikinks are able to propagate along three 
branches of the cactus meeting at the original defect tri- 
angle. Fig. 4(d). Thus a defect triangle is "made" from 
two antikinks with total spin 0. A state of two antikinks 
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FIG. 4: A state with a single defect triangle (a) evolves into 
a state with one longer-range bond (b), which can be viewed 
as two antikinks, shown as open and filled blue dots (c), with 
a total spin 0. State (d) is obtained from state (c) by an 
adiabatic exchange of the antikinks. Shaded triangles mark 
allowed paths of the antikinks. 



with a total spin 1 is also of interest to us, so we will con- 
sider a generic case of two antikinks with some specified 
spin projections, Fig. 4(c). 

The Husimi cactus is a toy model in which spinous 
have just enough freedom to exchange their locations by 
hopping along the three branches. That allows us to 
determine their quantum statistics by examining the re- 
sulting Berry phase acquired by the spinous. The two- 
spinon wavcfunction turns out to be antisymmetric under 
exchange, so antikinks are fcrmions. 

We first sketch an informal argument along the lines 
of Arovas et al. [21]. Starting with the state depicted in 
Fig. 4(c), we perform an adiabatic exchange of the two 
antikinks using the pivoting rule. We arrive at the state 
shown in Fig. 4(d) with dimers in the same positions but 
with one dimer reversed. If we use the same dimer basis 
to describe both states, the exchanged state acquires an 
extra factor of —1. The extra factor appears for any 
initial positions of the antikinks. 

The informal argument is supported by an explicit de- 
termination of the ground state of two antikinks on the 
cactus [18]. We find a nondegenerate, symmetric spatial 
wavefunction for total spin S" = and a doubly degen- 
erate, antisymmetric spatial wavefunction for S = 1, as 
one expects for two fermions. The ground state in the 
S = \ sector has energy J/2, twice the minimum energy 
of a free antikink. In contrast, two antikinks with 5 = 
form a bound state whose energy lies 0.06J below the 
bottom of the two-spinon continuum. The bound state 



has a small size f « 2.8 lattice units of the A chain. 

Curiously, antikinks in the 5 = ground state are 
bound not only to each other, but also to the original 
defect triangle: a pair is localized. The localization is 
topological in origin: like on the A chain, two antikinks 
cannot move through the same branch of the cactus: an 
antikink may only be followed by a kink. 

Let us discuss implications for the Heisenberg model 
on kagome. The cactus can be viewed as a tree made of 
triangles. Extrapolating results from a tree to a periodic 
lattice (even with the same coordination number) is only 
warranted for those physical quantities that are not sen- 
sitive to the long-distance properties of the lattice. For 
example, the band structure of a tree is quite different 
from that of a periodic lattice. Fortunately, the small 
size of a spinon pair means that its internal structure is 
determined by the local geometry of the lattice, so that 
the pair binding energy is expected to be roughly the 
same on the cactus, kagome, and hyperkagome. 

Both on kagome and hyperkagome one in four trian- 
gles carries a spinon pair. We are now facing a many- 
body problem of fermionic spinous with a concentration 
of 1/3 per site. The existence of spinous with Fermi-Dirac 
statistics in this and other frustrated two-dimensional 
magnets was conjectured previously [11, 22, 23] on the 
basis of a large- TV generalization of the Heisenberg model 
[24]. In these theories, spinous interact with an emergent 
lattice compact U(l) gauge field [2")]. We exploited an 
arrow representation of dimer coverings on kagome [C, 7] 
to define a fictitious U(l) gauge field living on links of a 
honeycomb lattice, Fig. 1(b). The arrows depict a quan- 
tized electric field of unit strength. Medial sites carry 
background charge Q = — 1, which is neutralized on av- 
erage by the charges of antikinks, Q = +2, and their 
pairs, Q = -1-4; kinks have charge Q = —2. It is note- 
worthy that kinks and antikinks are indeed expected to 
carry a small real electric charge [20]. 

The gauge-theory connection sheds light on the ori- 
gin of one-dimensional trajectories of spinous on kagome. 
An object of charge Q passing through a link increments 
the electric flux in it by Q, so antikinks can only move 
in the direction of electric flux. Fig. 1(b). As a result, 
spinous are constrained to move along one-dimensional 
paths that may be infinite (as on the Husimi cactus) or 
may terminate in a loop. A finite length of an antikink 
trail R would raise the triplet excitation energy of a pair 
by Ec ~ TT^ J/2i?^ [18]. Since R depends on the arrange- 
ment of dimers and may even vary from one spinon pair 
to the next; we only give an upper bound for the confine- 
ment energy. On kagome, i? > 10 so that < 0.05J. 
Thus the spin gap A = Ei, + Ec may vary from 0.06J to 
as much as O.llJ, depending on the valence-bond pattern 
in the ground state. 

Localization of antikink pairs can be traced to their 
interaction with the gauge field, whose electric flux is 
constrained to take on values ±1. Two antikinks travel- 
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ing along the same line would leave behind links with an 
unphysical electric flux ±3. Therefore, in order to move, 
the fermions making a pair must follow different paths. 
The tunneling amplitude is suppressed by a factor of or- 
der e^^/^, where ^ ~ 2.8 is the pair size and L is the 
length of the shortest loop on the medial lattice: L = 6 
for kagome, 10 for hyperkagome, and oo for the Husimi 
cactus. (The smallness of e~^/^ also makes it possible 
to treat heavy spinon pairs as stationary objects on time 
scales relevant to the motion of light spinous.) Singh 
and Huse [10] found that the presence of other pairs in- 
creases the amplitude of pair tunneling. Thus pairs form 
small localized clusters resonating around hexagons and 
preserving valence-bond order elsewhere. 

While our work agrees with the large-A^ approaches 
[11, 22, 23, 24] on the big picture of fermionic spinous 
interacting with a U(l) gauge field, there are important 
differences. First, the spinous are not free: exchange- 
mediated attraction binds them into small bosonic pairs. 
Second, the U(l) gauge field manifests itself as a quan- 
tized electric field, rather than a background magnetic 
flux, and frustrates the motion of spinon pairs. The mo- 
tion of spinous, which tends to scramble valence-bond 
order, is strongly suppressed as a result. These observa- 
tions lend support to the picture of a valence-bond crystal 
proposed by Marston and Zeng [9, 10, 11]. 

Aside from inelastic neutron scattering, another way 
to break spinon pairs is to apply a magnetic field strong 
enough to fully polarize antikink spins, so that the mag- 
netization reaches 1/3 of the maximum value. The liber- 
ated fermions will move at high speeds, inducing strong 
fluctuations of the electric field and restoring transla- 
tional invariance. A mean-field theory of spinous moving 
in the background of the emergent magnetic flux may be 
a good starting point for this partially polarized spin liq- 
uid. That would be interesting in light of recent numer- 
ical evidence for an incompressible quantum spin liquid 
at 1/3 of full magnetization [27, 28]. 
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